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We study the possibility for generating a new type of spin-orbit coupling for the center-of-mass 
motion of cold atoms, using laser beams that resonantly couple N atomic internal ground states 
to an extra state. After a general analysis of the scheme, we concentrate on the tetrapod setup 
(N = 4) where the atomic state can be described by a three-component spinor, evolving under 
the action of a Rashba-Dresselhaus-type spin-orbit coupling for a spin 1 particle. We illustrate a 
consequence of this coupling by studying the negative refraction of atoms at a potential step, and 
show that the amplitude of the refracted beam is significantly increased in comparison to the known 
case of spin 1/2 Rashba-Dresselhaus coupling. Finally we explore a possible implementation of this 
tetrapod setup, using stimulated Raman couplings between Zeeman sublevels of the ground state 
of alkali-metal atoms. 

PACS numbers: 37.10.Vz, 37.10.Jk 



I. INTRODUCTION 

The electron's spin degree of freedom plays a key role in 
the emerging area of semiconductor spintronics [l-Q . A 
first scheme for a semiconductor device is the spin field- 
effect Datta-Das transistor (DDT). It was proposed 20 
years ago Q and implemented recently Atomic and 
polaritonic analogs of the electron spin transistor have 
also been suggested An important ingredient of 

the DDT is the spin-orbit coupling of the Rashba [8l4l0| 
or Dresselhaus [ill) E3 types. This Rashba-Dresselhaus 
(RD) coupling scheme is described by a vector potential 
which can be made proportional to the spin-1/2 operator 
of a particle within a plane [13]. It applies to electrons 
^MjS E3] or atoms with two relevant internal states [lS - 

In the case of atoms, the spin-orbit coupling can be 
generated using two counterpropagating light beams [l7i — 
[H, UK (or two standing waves 0, HflUll) and a third 
beam propagating in an orthogonal direction, the beams 
being coupled to the atoms in a tripod scheme [l5l . l23l . l24j . 
The tripod atoms have two degenerate internal dressed 
states known as dark states, which are immune to atom- 
light coupling. The center-of-mass motion of the dark- 
state atoms is described by a two-component spinor and 
is equivalent to the motion of a spin-1/2 particle with 
spin-orbit coupling [i"6rl22| of the RD type. 

In the present article we investigate the possibility to 
generalize the RD spin-orbit coupling scheme to spins 
larger than 1/2. We show that this can be achieved us- 
ing cold atoms with more than two internal dark states. 
We start our analysis with the general scheme in which 
N laser beams couple N atomic internal ground states to 
a common excited state, thus forming the A^-pod setup 
shown in Fig.[TJ In the (AT+l)-dimensional Hilbert space, 
we identify N — 1 dark states, that is, zero-energy eigen- 




FIG. 1: iV-pod configuration. An atomic state |0) is coupled 
to N different atomic states \j), j = 1, N by N resonant 
laser fields. 



states of the atom-light Hamiltonian that are superposi- 
tions of the N ground states and are immune to atom- 
light coupling. 

Subsequently, we analyze the tetrapod case (N = 4) for 
which the center-of-mass motion of the dark-state atoms 
is described by a three-component spinor and thus corre- 
sponds to the motion of a spin-1 particle. We show that 
the resulting spin-orbit coupling can be made of the RD 
type and yields three cylindrically symmetric dispersion 
branches. Two of them are similar to those for the famil- 
iar RD spin-1/2 Hamiltonian, so the atom can exhibit the 
well-known quasirelativistic behavior [l8], [2(| for small 
wave vectors. Furthermore there is an extra branch with 
a flat dispersion around zero momentum. The formation 
of the latter branch leads to interesting phenomena, such 
as a possibility to have a negative refraction at a potential 
step, characterized by a larger amplitude as compared to 
the spin-1/2 case. 

Finally we explore a possible implementation of the 
tetrapod scheme for alkali-metal atoms using Raman 
transitions. To avoid a strong heating due to spon- 
taneous emission, all the states forming the tetrapod 
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scheme are chosen among the Zeeman sublevels of the 
atomic ground state, and are coupled by far-detuned Ra- 
man lasers beams. 



II. THE TV-POD SCHEME 
A. Atomic Hamiltonian 

We are interested in the center-of-mass motion of 
atoms in the field of several light beams. The atoms 
are characterized by N internal ground states |1), |2), ... 
, | AT), which are resonantly coupled to an extra state |0) 
by laser beams. This provides the TV-pod configuration 
shown in Fig. [TJ Note that the state |0) does not nec- 
essarily represent an electronic excited level; it can be a 
sublevel of the atomic ground state coupled to the states 
|1), 1 2), ... ,\N) via stimulated Raman transitions. A 
more detailed discussion on practical implementation is 
presented in the Sec. IVIIII 

The Hamiltonian describing the motion of an atom in 
the presence of the light beams is 



Ho = £-+V + Vi 
2m 



(1) 



where m is the atomic mass and p = —iftSI the atomic 
momentum operator. The terms Vq and V\ describe 
the atom-light interaction in the A^-pod configuration 
and a possible additional external potential, respectively. 
We assume for simplicity that all couplings \j) f-> |0), 
j = 1,...,N are resonant, so that Vb reads using the 
interaction representation and the rotating wave approx- 
imation: 



A' 



F =/!^%(r)|0><j|+H.c. 



(2) 



where £lj is the Rabi frequency that couples the internal 



state \j) to the common state |0), with j = 1,2, 
The coupling Vq can be rewritten as 



with 



V =m(r) \0)(B(r)\ + \B(v))(0\ 



1 N N 



(3) 



(4) 



Here \B) is the so-called bright (coupled) state and is 
the total Rabi frequency. 

The diagonalization of the atom-light interaction po- 
tential Vo is straightforward: 

(a) The coupling between the bright state \B) and the 
state |0) with a strength equal to the Rabi frequency f2 
in Eq. ^ gives rise to the two eigenstates 



± 



(|S>±|0»/V2, 



with energies ±M7. 

(b) The remaining orthogonal (N — l)-dimensional 
subspace corresponds to dark states. We denote \D n ), 
n = 1, . . . , N — 1 an orthonormal basis of this sub- 
space. All dark states are eigenstates of the Hamilto- 
nian Hq with zero eigenenergy: Hq\D h ) = . They 
are orthogonal to the bright state and to the state |0): 
(B\D n ) = (0\D n ) = 0. 

Although the eigenenergies of the dark states are 
position-independent, the states \D n ) depend on the 
atomic position through the spatial variation of the Rabi 
frequencies Q j . This leads to the appearance of the gauge 
potentials to be considered next. 



B. Adiabatic motion of dark-state atoms 

We now suppose that the atoms are prepared in the 
dark-state subspace, and that they move sufficiently 
slowly to remain in this manifold. This adiabatic approx- 
imation is justified if the light fields are strong enough, 
so that the energy difference ± hil between the dark-state 
manifold and the other eigenstates |±) of Vq is large com- 
pared to the detuning due to Doppler shifts. The atomic 
state- vector |$) can then be expanded on the dark-state 
basis 



JV— l 



(6) 



where \&j(r) is the wave function for the center-of- 
mass motion of the atom in the jth dark state. The 
atomic center-of-mass motion is described by an (AT — 1)- 
component wave function 



N-l 



■ , N. obeying the Schrodinger equation 



■* d T 



with the Hamiltonian 



H 



—{-mv 



A) 2 + $ + ^. 



(7) 



(8) 



(9) 



(5) 



The potentials governing the atomic center-of-mass mo- 
tion A, $, and V are (N — 1) x (N — 1) matrices. Here 
A and $ are the geometric potentials that emerge due to 
the spatial dependence of the atomic dark states [l^, [25l — 
[29j . The matrix A(r) represents a non-Abelian vector 
potential, with the matrix elements 

A n . m = ih(D n (r)\VD m (r)) , n, m = 1, . . . , N - 1 . 

(10) 

The matrix $(r) is an effective scalar potential known as 
the Born-Huang potential. It can be expressed through 
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the matrix elements of the vector potential between the 
dark states and the bright state \B) = \D ): 

$n, m = -^—A nfi • A , m , n, m = 1, . . . , N - 1 . (11) 
Am 

The matrix V(r) represents the restriction of Vi(r) to 
the dark state subspace. For simplicity we assume in the 
following that (i) the matrix elements of Vi between the 
dark-state manifold and the states \B) or |0) are negli- 
gible, so that V\ cannot cause any significant departure 
of atoms from the dark-state manifold; (ii) V is propor- 
tional to the identity matrix in the dark state subspace, 
so that it does not break the gauge symmetry of (A, $). 
For the particular case of alkali-metal atoms, this oc- 
curs when the trapping is provided by far-detuned laser 
beams. The confinement potential is then the same for 
all sub-levels of the electronic ground state, in particular 
for the states \j) (j = 1, . . . , N) considered here. 

The non-Abelian vector potential A provides a curva- 
ture (or effective "magnetic" field) 

B = VxA + 4rAxA. (12) 
in 

The first term represents the usual curl. Note that the 
second term A x A does not vanish in general, since 
the Cartesian components of the vector potential A do 
not necessarily commute (i.e. the vector potential is non- 
Abelian) . Therefore in contrast to the Abelian case, even 
a constant vector potential can produce a nonzero curva- 
ture and thus provide nontrivial topological effects, lead- 
ing, for example, to unusual dispersion curves. 

III. EFFECTIVE FIELDS GENERATED BY 
PLANE- WAVE LASER BEAMS 

A. Dark states and gauge potentials 

From now on we focus on the case where the laser 
beams represent plane running waves characterized by 
wave vectors k,, j = 1, . . . , N. We suppose that the N 
Rabi frequencies have equal amplitudes and read 

a = ^Li2e^ r , j = 1,2,...,A. (13) 

At this stage the directions of the wave vectors kj are 
still arbitrary; we will address some specific geometries 
in Sees. |V] and ED 

A convenient orthogonal set of N — 1 normalized dark 
states is 

1 N 

= 7 =5>V 2 " WA, - !;krr , (14) 
V A j=1 

with n = 1,2,..., A— 1. Note that the bright state given 
by Eqs. gj) and ([T3J) has the form of Eq. (TTJ]) with n = 0, 
so we will use in the following the notation \Dq) = \B). 



Equations (|10l) and (p~4|) provide the matrix elements 
of the vector potential 

h N 

It is evident that the vector potential A„ jm depends only 
on the difference n — m, i.e. A„ )OT = A„_ m 

Since the vector potential given by Eq. (Tl5)) is constant 
in space, the effective magnetic field (IT21) simplifies to 
ihB = A x A. Using Eq. (fl~5|) . it can expressed in terms 
of the off-diagonal matrix elements of the vector potential 
A„ ; o and A 0)TO : 

B n , m = -rAn.o x Ao, TO . (16) 



B. Vector potential and angular momentum 

We now address the following question: Can the vector 
potential A be made proportional to a three-dimensional 
(3D) angular momentum operator J, that is, A = 7J , 
where 7 is a constant? If the answer was positive, this 
would allow one to achieve a three-dimensional RD-type 
coupling. This would be formally similar to the effective 
spin-orbit interaction discussed in [30f , arising from non- 
Abelian gauge fields in molecular physics. However as we 
see now, one cannot use the present scheme to achieve 
AaJ. 

The angular momentum operator is known to obey the 
following relations: 

J x J = ihJ . (17) 

If A = 7 J , the cross product of the vector potential 
should be proportional to the vector potential itself: 
7 A x A = ihA or simply 7B = A. Using Eq. (fl"6j) . 
the last relationship would lead to 

7A„ i0 x A^ q = -ihA n - mfi , n, m = 1, . . . , N - 1 . 

(18) 

Multiplying Eq. ([18]) by A* rn , the left-hand side of the 
resultant equation is zero. Thus one arrives at 

A^ • A„_ m , = 0. (19) 

Equation (Ti"T)l) should hold for all possible values of n and 
m. In particular, by taking m = 1 and n — 2, one finds 
Al • A^o = 0. This equation can be fulfilled only if 
A^o = 0. Then by taking m = 1, the relationship (TT8"|) 
yields that A„_i.o = A n+PtP+ i — for integer n and p. 
This means that the vector potential A = 7J should be 
identically equal to zero. 

In this way, we have proved that when using the N- 
pod scheme with plane waves of equal amplitudes it is not 
possible to generate a nonzero vector potential which is 
proportional to the 3D angular momentum operator J. 
In other words, it is not possible to produce a 3D spin- 
orbit coupling of the RD type using the A-pod scheme. 
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Yet one can get a two-dimensional (2D) RD coupling by 
means of the N-pod scheme. This includes not only the 
usual spin- 1/2 RD coupling but also a generalized 2D RD 
coupling for the spin-1 case, as we shall see later on. 



IV. PLANE MATTER- WAVE SOLUTIONS 

We suppose in the following that the external poten- 
tial V is uniform in space. In this case the Schrodinger 
equation ([8j has plane-wave solutions: 



$kM) = * k e 



ik-r— uj-irt 



(20) 



where the dimensionless parameter /3 = /3k depends gen- 
erally on the wave- vector k. The corresponding eigenval- 
ues of the Hamiltonian are 







2m 



(k + /3k) 2 + (2 - (3 2 )k< 



V. 



(27) 



For k = all the eigenenergies fkv^ are equal and do 
not depend on the branch parameter j3. Consequently all 
dispersion branches merge to luq = ujq at the origin where 
k = 0. To find the eigenstates and the eigenenergies for 
k 0, one needs to specify the arrangement of the wave 
vectors k,. 



where Wk is an eigenfrequency and 'I'k is a k-dependent 
spinor: 



(21) 



(AT-l),k 



Note that the direction of the wave vector k is arbitrary 
and it is not related to the wave vectors of the light beams 
ki. 

The k-dependent spinor 'I'k obeys the stationary 
Schrodinger equation 

i? k * k = fta; k *k, 
with the k-dependent Hamiltonian 



/4 = _ A A ■ k + ^a 2 + $ + v 



2m m 



2m 



(22) 



Exploiting Eqs. (ITU and (fT5|). the scalar term A 2 /2m+<I> 
takes the form 



A 2 + $ = > k . 

2m ) n rn 2m N ^ > 



(23) 



If the wave vectors of all the Rabi frequencies have the 
same modulus k 2 = 2k 2 , the term 



1 » 2 * h2 K 2 t 
- — A + $ = / 

2m m 



(24) 



is proportional to the unit matrix / for any arrangement 
of the wave vectors (both planar and 3D). In this case 
the Hamiltonian (f2"2"j) simplifies to 



iTk = tt~ I Kk — 2A • k + 2hn z )+V. 
2m 



(25) 



If the external trapping potential V is proportional to 
the unit matrix, the eigenvectors of the Hamiltonian 
Hk are also the eigenvectors of the operator = A-k/fc 
representing the projection of the vector potential along 
the wave vector, 



(26) 



PLANAR GEOMETRY 



A. Wave vectors on a regular polygon 

Let us analyze a situation where the wave vectors k, 
are situated in a plane and form a regular polygon 



kj = \/2/c[— cos a.je x + sin ay e^] 



'e_ 



(28) 
(29) 



with e± = -^( e x i i e y), where ay — 2nj /N is the angle 
between the wave vector and the x axis. The scalar and 
vector potentials, Eqs. (fTTj) and (|T5|). take the form 



2^.2 



2m 



■(Sm,lS n ,l + $m,N-l8n,N-l) j (30) 

(31) 



-/IK ^ e-fc^n,m±l 



The vector potential is thus a tridiagonal matrix whose 
elements are proportional to e x ± ie y , whereas the scalar 
potential $ n ,m is a diagonal matrix with nonzero ele- 
ments only for n = m = 1 or n = m = N — 1. 

Note that the matrices A x and A y , are proportional 
to the x and y components of the angular momentum 
operator J only for the tripod (N — 3) and tetrapod 
(N = 4) schemes. In these cases the scalar potential is 
proportional to J 2 . 

The projection of A„^ m along the wave vector is 

(A k )„, m = --^[6 n , m+ ie^ + ^ vm _ie- 4V J , (32) 

where (p is the angle between the wave vector k and the 
x axis. The eigenvectors of this operator are 



ZL£) 



sm 



sin (2f ) e«*> 



V sin ((/V- 1)^2) e l ( w -2)^ 



(33) 



with g = 1, ... N — 1. The corresponding eigenvalues are 
given by Eq. (|26|) with 



(34) 
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Q 2 ~ e ifc " 



FIG. 2: (Color online) Planar arrangement of laser beams for 
tripod (a) and tetrapod (b) setups. 



It is to be emphasized that the dimensionless parameter 
/3 does not depend on k for this particular geometry. 
The vectors "J^ represent eigenstates of the Hamiltonian 
with eigenenergies given by Eqs. (|27p and (|34l) . This 
provides N — 1 dispersion branches. 



B. Tripod setup 

Consider first the tripod setup (N — 3) in which the 
wave vectors kj form an equilateral triangle [Fig. [UJa)]. 
The parameter /3 featured in Eqs. (J25]) and (jM)) then 
takes the values h(3/y/2 = ±ft/2, representing the eigen- 
values of the projection of a spin- 1/2 on a given axis. In 
such a situation the operator A is related to the spin 1/2 
operator hcr±, providing the RD coupling along the xy 
plane as in the previous studies |16l - l22j |: 



A = — finer i 



(35) 



It is noteworthy that the present setup produces a cylin- 
drically symmetric spin-orbit coupling in a more straight- 
forward manner than the previously suggested tripod 
schemes. Those schemes involved two counterpropa- 
gating light beams [i~7rfi~9l [22[ (or two standing waves 
[la . [20, [21[) and a third beam propagating in an orthog- 
onal direction. Consequently, one needed to add a detun- 
ing potential and make the amplitudes of the Rabi fre- 
quencies asymmetric in order to have dispersion curves 
of the RD-type, with the proper cylindrical symmetry 
[l6l - |22l | . On the contrary, for the present regular polygon 
arrangement of wave vectors, the dispersion relation is 
naturally symmetric as long as the amplitudes of all four 
Rabi frequencies are equal. 



C. Tetrapod setup 

For N — 4 one arrives at the tetrapod setup involv- 
ing two pairs of counterpropagating laser fields shown in 
Fig. [H[b). In this case the vector potential reads 








-e. 



ie.„ 



(36) 



C0 k /(0 




FIG. 3: (Color online) Dispersion curves for the tetrapod 
scheme calculated using Eq. (|40p for V = 0. Here uq = 
Fin 2 /m. 



The possible values for the parameter featured in 
Eq. (|26| are h(3 — 0, ±h, representing the eigenvalues 
of the component of a spin 1 along a given axis. Conse- 
quently the operator A is proportional to the projection 
J_l of a spin 1 operator along the xy plane 



Jx^x 



J y e y . 



(37) 



The scalar potential can be represented in terms of the z 
component of the spin operator 



h 2 n 2 
2m 



-J 



(38) 



The eigenstates and the eigenenergies of the Hamilto- 
nian are now 




k V2 



-1 


J2ia 



(39) 



and 

fiwf = - — yk 2 + 2nk(3 + 2k 2 j + V , /3 = 0,±1. (40) 

For /3 = ±1 the dispersion curves shown in Fig. [3] are 
analogous to those of the spin-1/2 RD model. An addi- 
tional dispersion curve with /3 = represents a parabola 
centered at fc = 0. 

The dispersion curve with (3 — —1 has its minimum 
at hjj — h 2 K 2 /2m, whereas the other two dispersion 
branches have minima at the double energy tuu = H 2 k 2 /m 
(for V — 0). Therefore, all dispersion curves have a 
strictly positive minimum energy. This nonzero mini- 
mum originates from the micromotion of the atom in the 
light field, caused by nonadiabatic transitions between 
the dark and bright states pJl], HH . The associated ki- 
netic energy gives rise to the scalar potential given by 
Eq. (I30|) . which has a nonzero contribution even when 
acting on the dark states. 
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Finally, we note an important difference in the "topol- 
ogy" of the eigenfunctions for the RD spin-1/2 and spin-1 
problems, even thought the f3 — ±1 branches have the 
same dispersion in the two cases: The wave functions "f^ 
exhibit a 7r Berry's phase in k space in the spin-1/2 case, 
whereas this Berry's phase is absent for the spin-1. 



Hence, one cannot generate a 3D Hamiltonian with RD- 
type spin-orbit coupling for all directions of the atomic 
motion. This is a particular case of the general conclusion 
reached in the Sec. IIII B[ 

VII. TRANSMISSION BY A POTENTIAL STEP 



VI. TETRAHEDRON GEOMETRY 

In this section we present an example of a nonplanar 
setup, which has some advantages with respect to the 
planar configuration investigated in the previous section, 
because it leads to a simpler scalar potential. We consider 
again the tetrapod setup (N — 4) with wave vectors kj 
arranged in a regular tetrahedron geometry: 



3 ^ f ■ 



(41) 



where k, = kj/fcj is a unit vector. More precisely, we 
choose 

k 1>3 = K'(±e y V2-e z ) , k 2A = K'(±e x V2 + e z ) . (42) 
Using Eq. (fl"5j) . the vector potential then reads 



V2 







V2e 



y/2e z 

: + 





(43) 



For atoms moving in the xy plane the vector potential 
can be expressed in terms of a spin-1 operator in the xy 
plane: = — k'3±. Hence we obtain as before a RD- 
type spin-orbit coupling for the atomic motion in the xy 
plane, characterized by the dispersion relation shown in 
Fig. |31 Yet we are now dealing with a 3D problem, so 
the same dispersion also characterizes the atomic motion 
along two other planes perpendicular to the vectors e x + 
e y and e x — e y . By making an atomic lattice along these 
directions, the atomic tunneling will be influenced by a 
spin-1 RD coupling, thus extending the previous studies 
of spin-1/2 RD coupling in lattices [33|. This will be 
investigated in a separate study. 

A distinguished feature of the tetrahedron geometry 
is that the scalar potential is proportional to the unit 
matrix /: 



ftV 2 f 
<t> = 1 

2m 



(44) 



Thus for atoms placed in a 3D lattice, there is no energy 
mismatch between different dark states located in adja- 
cent sites. This contrasts with the planar tetrapod case 
(Eq. ([55])). where the spin components are likely to get 
frozen in the lattice because tunneling matrix elements 
are normally much smaller than the atomic recoil energy, 
which gives the scale for the scalar potential. 

It is noteworthy that the z component of the vector 
potential given by Eq. (|4"3")l is not proportional to J z . 



A spectacular consequence of spin-orbit RD coupling is 
the negative refraction and reflection that occurs when a 
matter wave is incident on a potential step. The problem 
was investigated for spin-1/2 atoms [H,|34| and electrons 
35]. In this case one can calculate relatively easily the 
transmission and reflection of the atomic wave packet. 
For small wave vectors of the incident atoms, k <C K, 
the transmission probability is close to unity at zero an- 
gle of incidence. Here the parameter k characterizes the 
strength of the spin-orbit interaction, see Eq. ([33]) . This 
nearly complete transmission is a manifestation of the 
Klein paradox appearing also for electron tunneling in 
graphene [36j . For a nonzero angle of incidence, the 
transmission probability is less than 1 and decreases with 
increasing angle. Furthermore the transmitted matter 
wave experiences negative refraction [l8| . similar to the 
case of electrons in graphene [37j ■ 

Particles with a spin larger than 1/2 have additional 
degrees of freedom, which modifies the continuity condi- 
tions at the potential step. This can lead to a significant 
increase of the transmission probability of atoms, as we 
show now for particles submitted to a spin-1 RD cou- 
pling. 



A. The Hamiltonian 

We consider in this section the motion of a particle in 
the xy plane described by the Hamiltonian 

H = — (p 2 + 2hnp ■ J ± + 2h 2 K 2 ) + V(x) (45) 
2m 

where = J x e x + J y e y is the projection of spin-1 
operator onto the xy plane. Such a Hamiltonian can 
be obtained using the tetrapod setups described in the 
Sees. IV CI and I VII The external potential V(x) is given 
by the step function along x 



V(x) 



0. 



x < 



Vo, x>0 



(46) 



with Vo > 0. It is convenient to introduce the wave vector 
fco = 2mVo/h 2 n characterizing the height of the barrier. 

For a constant potential the eigenvalue equation has 
plane-wave solutions (|20[) characterized by the spinor 
part ^ [Eq. (|3"5j) ]. The corresponding eigenvalues Hco^. 



are given by Eq. (|4"01) with k — ^jk x + fc 2 and are plot- 
ted in Fig. [3J Additionally there can be evanescent wave 
solutions localized in the vicinity of the potential step in 
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FIG. 4: (Color online) Wave numbers of reflected and trans- 
mitted waves and energy conservation at the potential step. 



&3 k 2 ki 




FIG. 5: (Color online) Reflection and transmission of atoms 
at a potential step. In addition there is an evanescent trans- 
mitted wave with kg = k v e v + iqe x . 



the x > region. In that case we have k x — iq with 
q > 0, giving 



k=Jk 2 -q 2 , q 2 <k 2 y . 



(47) 



For the present problem, only the evanescent wave with 
(3 = will play a role, 



k V>1 k y>1 \ q+k v 

q-k y 



(48) 



where ct is the normalization factor. 

r»y ,q 



B. Incident waves with j3 = I 

In this paragraph we restrict our analysis to the case 
where the incident atom is prepared in the upper disper- 
sion branch (/? = 1) in the region x < 0. Denoting its 
wave vector by k, the incident wave is 



*in = *£e ikr 



(49) 



The potential step is assumed to be high enough, 

k Q >k 2 /K + 2k, (50) 



so that there can be no propagating transmitted waves 
with chirality /3 = or j3 = 1 (see Fig. 2]) . At the same 
time, to allow for propagation of plane waves in the region 
x > for the lower dispersion branch (3 = — 1, the step 
height should not be too large: 



fco < k + k 2 1 n + 2k . 



(51) 



The directions of reflected and transmitted waves are de- 
picted in Fig. [5] The reflected waves generally contain 
all three components, 

*rcfi = ri^e* 1 " + r 2 *° 2 e lk2 r + r 3 * k V k3 r , (52) 



where k\ = k, k2 = \/k 2 + 2nk and k% = k + 2k. The 
reflection angles are n — a, ir — a 2 , and ir — a^, with 
a2 = arcsin[sin(a)fc/fc 2 ] and 0J3 = arcsin[sin(a)fc/fc3]. 
The transmitted waves are 



•t 6 *L e lker . 



(53) 



— fc K, k 5 = K + 

k 2 + 2nk — koK (see 



where ki — K — y/(k + n) 2 
y/ (k + k) 2 — koK , and fee = 
Fig. @|. The first and second transmitted waves ex- 
perience negative and positive refraction, respectively, 
and propagate at the angles ir — a 4 and CV5, where 
a 4 = arcsin[sin(a)fc/fc 4 ] and a 5 = arcsin[sin(a)fc/fc 5 ]. 
On the other hand, due to the condition ([501 the third 
transmitted wave with the helicity /3 = is an evanes- 
cent one along the x axis and thus is characterized by 



the wave vector k« = k v e v 



iqe x , with k y 



ks'ma and 



q 



h 2 



k 2 . Note that there is no evanescent trans- 



mitted wave in the upper dispersion branch (/3 = 1) be- 
cause it cannot comply with the momentum conservation 
along the interface in addition to the energy conservation. 

The multicomponent wave function and its first deriva- 
tive in the x direction are required to be continuous at 
the barrier (x — 0), providing six equations containing 
six unknown coefficients r l5 r 2 , r 3 , t 4 , t 5 , and t 6 . Of 
special interest is the situation where fco = 4fc. In this 
case the wave number of the first refracted wave coincides 
with the wave number of the incident wave, fc 4 = fc, so 
the angle of refraction is equal to the angle of incidence 
for the first reflected wave. a 4 = a. 

The analytical solution for the six coefficients is gen- 
erally complicated. It is instructive to obtain approxi- 
mate solutions for small wave vectors and small angles 
of incidence, fc <C k and a -C 1. In such a case one 
can restrict to reflected (|52|) and transmitted (|53l) waves 
containing only the contributions of ki, k 2 , k 4 , and k^. 
The transmitted wave with kg represents a rapidly decay- 
ing evanescent wave characterized by a spinor component 
given by Eq. (|48|) with q > k y : 



(54) 




The continuity of the wave function at x = gives 



„ik 2 r 
k 2 e 



(55) 
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In addition, we require continuity of the derivative in the 
x direction for the component with (3 = 0, which is the 
largest: 



k 2 cos(tt - a 2 )r 2 *k 2 = iqa^l 



k 6 



(56) 



k and a 2 ~ y / k/2nsin(a), with k 2 = 



Here k\ « fco/2 
\j2xk and q « ^/ k(/cq 



2k). Using the spinors (|39|) and 
54| one obtains the following solution to Eqs. (|55|) and 
5F 



*4 



'"2 



2 cos a 



i(a+a 4 ) 



cos a + cos «4 
cos «4 — cos a i2 Q 
cos a + cos «4 
Vfcp - 2k 
Vk + i^ko/2 - k 



cos a tan 



■ Ct4 



(57) 
(58) 

"(59) 



The calculated reflection and transmission coefficients r\ 
and ti obey the probability conservation up to terms of 
the order of 0(a 2 ): 



l9 C0S«4, .<- 

N 2 + 

cos a 



f 



(60) 



If the barrier height is such that ko = 4k, we have ay = a. 
In that case |i4| « 1 and |ri| « leading to an almost 
perfect negative refraction, at the exactl opposite refrac- 
tion angle, provided k <C k and the angles of incidence 
are not too large. 

Figure [5] presents the comparison of the transmission 
probabilities for the spin-1 and spin- 1/2 RD coupling us- 
ing the exact numerical solutions of the continuity equa- 
tions at the boundary x — 0. For the spin- 1/2 case the 
incident wave is also prepared in the upper dispersion 
branch. The figure shows a marked increase in the trans- 
mission probability for small angles of incidence in the 
case of spin 1. Note that the transmitted waves expe- 
rience negative refraction both for the spin-1 and the 
spin- 1/2 cases. 



VIII. IMPLEMENTATION OF THE TETRAPOD 
SETUP WITH ALKALI-METAL ATOMS 

We now discuss a possible implementation of the tetra- 
pod scheme. We consider the case of alkali-metal atoms, 
which are the most frequently used in current experi- 
ments. In order to avoid a strong heating due to sponta- 
neous emission, we study the case where the state |0) is 
actually one of the Zeeman sublevels of the ground state. 
The states \j) (with j = 1, ... ,4) are also Zeeman sub- 
levels of the ground state, and the coupling between the 
state |0) and a state \j) is provided by a pair of laser 
beams that induce a Raman transition under the con- 
dition of the two-photon resonance. The use of Raman 
transitions in this context is an extension to the tetra- 
pod case of a recent proposal 38] to implement a A-type 



T(a) 



1 

0.8 
0.6 
0.4 
0.2 


10 20 30 40 50 60 70 80 90 
a 

FIG. 6: (Color online) Transmission probability of negatively 
refracted atoms at a potential step as a function of the angle 
of incidence a for a spin-1 (solid red) and spin-1/2 (dashed 
green) systems. The parameters used for the calculation are 
k/n — 0.1 and ko/K = 0.4 for both systems. 




>0 



FIG. 7: Directions of the laser beams used for implementing 
the tetrapod coupling scheme with alkali-metal atoms, via 
stimulated Raman transitions. 



scheme for the generation of an effective magnetic field 
by means of the counterpropagating laser beams [lU, [3t| . 

We recall that the electronic ground level nS\j 2 of 
alkali-meal atoms is split by hyperfine interaction in 
two sublevels with angular momenta F = 1+1/2 and 
F = I — 1/2, where I is the nuclear spin. We consider 
in the following the case 1 = 3/2 that is relevant for 
lithium ( 7 Li, n = 2), sodium ( 23 Na, n = 3) or rubid- 
ium ( 87 Rb, n = 5). In order to minimize the rate of 
spontaneous emission processes, we restrict to Raman 
transitions that are far detuned from the resonance with 
the "true" excited states nPi/ 2 or nP 3 / 2 of the D\ or 
D 2 transitions. More precisely the typical one-photon 
detuning of the beams involved in the Raman process is 
chosen much larger than the hyperfine structure of the 
excited level nP 1 / 2 or nP 3 / 2 (0.8 GHz for the hyperfine 
splitting of the level 5Pi/ 2 of 87 Rb). At the same time 
the one-photon detuning should be smaller than the fine 
structure splitting, that is, the difference between the en- 
ergies of nP 1/2 and nP 3/2 (7000 GHz for 87 Rb). When 
the one-photon detuning exceeds the hyperfine splitting, 
the nucleus angular momentum does not play any role in 
the selection rules that determine the allowed transitions 
for photon absorption or emission. For the D\ (D 2 ) tran- 
sition, the allowed couplings are the same as between a 
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FIG. 8: Implementation of the tetrapod scheme for alkali- 
metal atoms with two hyperfine levels of angular momen- 
tum F = 1 and F = 2. The laser couplings involved in 
this scheme correspond to stimulated Raman transitions be- 
tween hyperfine states of the ground atomic level. We choose 
|0> = \F = l,m F = 0). (a) The laser beams A, A 1 , B and 
B' induce the transitions |0) — > |1) = \F = 2,m,F = 1} and 
|0) -¥ |3) = \F = 2,m F = -1). (b) The laser beams A" and 
B" induce the transitions |0) — » \2) = \F = 1,ttif = 1) and 
|0> ->• |4> = \F = l,m F = -1). 



spin-1/2 ground level and a spin-1/2 (3/2) excited level. 
In particular, the only allowed Raman transitions cor- 
respond to a change Amj = or Amj = ±1 of the 
azimuthal quantum number uij. 

A scheme that fulfills the aforementioned constraints 
is represented in Figs. [7] and [51 The atomic motion along 
the z direction is supposed to be frozen thanks to a trap- 
ping potential muj^z 2 /2 such that Hu) z is much larger 
than the atomic kinetic energy. The atom is placed in a 
uniform magnetic field Bq directed along the x direction. 
The role of this magnetic field is to allow for a selective 
Raman excitation between two given Zeeman sublevels. 
More precisely the Larmor frequency = fJ-BBo/h (/^b 
is the Bohr magneton) is chosen much larger than the 
two-photon (Raman) Rabi frequency f2. Typically we 
choose wl/27t on the order of a few MHz (i.e. Bo on the 
order of a few Gauss) and fi/27r in the range 10 5 — 10 6 Hz. 
The latter choice is sufficient to ensure that the split- 
ting Ml between the dark-state manifold and the states 
|±) is large compared to the two-photon Doppler shift, 
as required for the adiabatic approximation to be valid. 
The state |0) is chosen equal to the \F — l,mp = 0) 
sublevel and the states \j) with j — 1,...,4 are the 
\F = 2, nip = ±1) and \F — 1,itlf = ±1) sublevels. Here 
the quantization axis is the x axis, parallel to the direc- 
tion of the magnetic field Bo. As indicated in Fig. [SJa), 
the transition between |0) and |1) = \F = 2, mp = +1) is 
driven by a pair of laser beams [A, B) with a frequency 
difference equal to Whf + where lom is the hyper- 

fine splitting between the F = 1 and F = 2 manifolds 
(whf/27r is on the order of 7 GHz for 87 Rb). The laser 
beam A propagates along the y axis (wave vector fce y , 
where e y is a unit vector). It is linearly polarized along x, 
so that it carries no angular momentum along the x axis. 
The laser beam B propagates along the x axis (wave vec- 
tor — fce x ) and is circularly (<7-) polarized. In the tran- 



sition |0) — > |1) the momentum change of the atom is 
tiki = hk(e x + e y ). One can readily check that the tran- 
sition |0) — > |1) is the only one that is driven resonantly 
by this pair of beams, thanks to the fact that the Lande 
factors are opposite for the F — 1 and F = 2 manifolds, 
as one can see in Fig. [5Ja). Similarly, the transition be- 
tween |0) and |3) = \F = 2, mp = —I) is driven by a 
pair of laser beams (A',B') with a frequency difference 
equal to Whf — <^l/2. The beam A' propagates along y 
with wave vector — ke y and is linearly polarized along x. 
The beam B' propagates along x with wave vector fce x 
and is circularly (cr+) polarized. The atomic momentum 
change in the transition |0) — » |3) is Kks = — Kkx- The 
difference in the frequencies of A and A' is chosen large 
enough so that no transition is driven with a significant 
probability by the pairs of beams (A, B') and (A' , B). 

The two remaining states of the tetrapod configura- 
tion are |2) = \F = 1, m F = +1) and |4) = \F = l,m F = 
— I). The coupling between these states and the state |0) 
is provided by a single pair of laser beams (A", B"), as in 
the recent experiment [38] , in which the A (ladder) type 
coupling was generated within the Zeeman sublevels of 
the F — 1 manifold. The wave vector of A" is fce y and 
this beam is linearly polarized along x. The beam B" 
propagates along x with wave vector fce x and is circu- 
larly (<r + ) polarized. The frequency difference between 
the beam A" and B" is w L /2 so that the pair (A",B") 
resonantly drives the transition |0) — > |2) with a mo- 
mentum transfer ftk.2 = hk(e x — e y ), and the transition 
|0) — » 1 4) with a momentum transfer ftk.4 — — fi,k2- Note 
that here again we take advantage of the different signs of 
the Lande factors of the F = 1 and F = 2 manifolds: The 
pair of beams (A" , B") cannot resonantly drive a tran- 
sition between two sublevels of the F = 2 manifold []see 
Fig. [FJb)]- Note also that another consequence of two- 
photon processes is a modification of the energies of the 
states via the absorption and stimulated emission of 
photons in the same laser beam. It can be accounted for 
by including these energy shifts in the choice of the two- 
photon detunings and, for example, by taking advantage 
of the (small) second-order Zeeman shift. 

This configuration therefore constitutes a suitable im- 
plementation of the scheme discussed in the first part of 
this article. The momentum transfers Mcj = fi,fc(±e x ± 
e y ) form a square in the xy plane shown in Fig. HJb) 
(subject to the rotation of the coordinate system by 45°). 
The intensities of the various beams can be adjusted so 
that all Rabi frequencies ilj are equal, once the Clebsch- 
Gordan coefficients associated to each Raman transition 
have been taken into account (note that the two-photon 
Rabi frequencies for |0) — > |2) and |0) — > |4) transitions 
are equal by construction). With a one- photon detun- 
ing of 3 nm , which represents 1/5 of the fine structure 
splitting for rubidium atoms, the residual photon scatter- 
ing rate is below ls^ 1 for a two-photon Rabi frequency 
H/(2tt) = 10 5 Hz. The corresponding heating rate is thus 
small enough to provide enough time for the investigation 
of the RD coupling studied in this article. 
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IX. CONCLUSIONS 

In this article we have explained how to produce a 
spin-orbit coupling of the RD type for a spin larger than 
1/2. Our scheme makes use of cold atoms with three or 
more internal dark states so that their quasi-spin is equal 
to or greater than unity. We have analyzed a general 
scheme in which N laser beams couple N atomic internal 
ground states to an extra state, thus forming an A^-pod 
setup of light-matter interaction. In this case the atoms 
have N — 1 dark states representing superpositions of 
the N ground states that are immune to the atom-light 
coupling. 

We have analyzed in detail the particular case of the 
tetrapod setup (N = 4), in which the center of mass mo- 
tion of the atoms in their dark state manifold is described 
by a three-component spinor and thus corresponds to 
the motion of a spin-1 particle. We have shown that 
the resulting spin-orbit coupling can be made of the RD 
type and yields three cylindrically symmetric dispersion 
branches. Two of them are similar to those known for 
the familiar RD spin- 1/2 Hamiltonian, so the atom can 
exhibit a quasirelativistic behavior [18- 20] for small wave 
vectors. Furthermore, we have shown that there exists 
an extra branch with a fiat dispersion around zero mo- 
mentum. We have studied the modifications that this 
extra branch brings to the problem of negative refraction 



of matter waves on a potential step, and shown that it 
enhances the negative refraction probability. 

Finally we have discussed a possible implementation 
of the tetrapod setup with cold alkali-metal atoms. We 
have shown that in order to avoid heating due to spon- 
taneous emission, it is possible to choose all the states 
involved in this tetrapod scheme among the various Zee- 
man sublevels of the ground atomic state. All laser cou- 
plings are then provided by stimulated Raman transi- 
tions. For rubidium atoms, realistic parameters yield a 
residual spontaneous emission rate below Is -1 , which 
makes the observation of this spin-orbit coupling scheme 
experimentally feasible. 
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